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We review the method developed in Pisa to determine the topological susceptibility in lattice QCD 
, and present a collection of new and old results obtained by the method. 
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1 Introduction 



This article is dedicated to the study of topological properties of QCD on the lattice. We 
will describe the method developed in Pisa to measure the topological susceptibility on 
the lattice which is usually known as field-theoretical method and we will review results 
obtained by the method together with some recent applications. 

Topology is a difficult subject to be studied on the lattice because from a strictly math- 
ematical point of view it cannot be defined on a discrete space. Its meaning is recovered 
in the continuum limit, where the physical scale gets well separated from the ultraviolet 
(U V) one. T he topological susceptibility enters in many contexts of QCD phenomenol- 

o 



°gy • ^ can De treated as usual expectation values of field operators in Quantum 

Field Theory. In this approach one does not try to assign a definite topological sector to 
each gauge field configuration. Instead appropriate renormalization constants are used to 
relate the lattice averages to the corresponding continuum ones. The method developed 
in Pisa consists in defining the renormalizations and in giving a prescription to compute 
them on the lattice, as one would do for any other field-theoretical expectation value. We 
will review this method in Section 2. 

In Section 3 we collect some results obtained in the past by using the field-theoretical 
^ ■ method in QCD at zero and finite temperature, with and without dynamical quarks. In 

Section 4 we present new results concerning the application of the method at finite tem- 
perature on anisotropic lattices and at finite density. A feature of these last results is that 
several determinations of the topological susceptibility are obtained with a fixed UV cut- 
off and this fact greatly reduces the number of renormalization constants required in the 
computation and makes the method even more practical and beautiful. 

We have written this article to celebrate the 70th birthday of Adriano Di Giacomo: 
the study of topology on the lattice is surely one of the leading interests in his research 
activity, an interest which both of the authors have been lucky enough to share with him. 

2 The topological susceptibility 

The topological susceptibility \ is defined as the zero momentum two-point function of 
the topological charge density operator q(x) 

X = = |d 4 x (q(x)q(O)) , Q = J d 4 x q(x) , (1) 
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q(x) = ^e^ pa F^(x)F^(x) . (2) 

It contains information about the dependence of the QCD free energy on the topological 
8 parameter around 8 = 0, 



X- T T J/12 



lnZ{6)\ g=Q , Z{6) = J[dA] e - /***Ajm>0»0 jOQ (3) 

and it enters in various aspects of QCD phenomenology by regulating the realization of 
the U(1)a axial symmetry. For instance, its value in the pure gauge theory is directly 
related, at the leading order in 1/N C (N c is the n umber of colours), to the mass of the r/ 

meson by the Witten-Veneziano mechanism^^l which predicts x ra (180 MeV) 4 . 

Topological properties of QCD are of non-perturbative nature: lattice simulations are 
therefore the natural tool to investigate them. On the lattice a discretized gauge invari- 
ant topological charge density operator qi{x) can be defined, and a related topological 
charge Qi = J2 x 1l{x), with the only requirement about the formal continuum limit 

lim a ^o Ql(x) = a 4 q(x) where a is the lattice spacing. A commonly used definition i s 6 

-1 ±4 

^( x ) = ^2 zZ W^ Tr ( n ^( x ) n p^( x )) i ( 4 ) 

where Yl fiu (x) is the plaquette operator in the \iv plane, e M ^ p(J is the Levi-Civita tensor 
for positive directions and is otherwise defined by the rule e^ pcr = —^i-fj) v pa- 

A proper renormalization must be performed when going towards the continuum limit. 
In spite of the formal limit, the discretized topological charge density renormalizes mul- 
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tiplicatively L '-' 2, (in presence of dynamical fermions also additive renormalizations can be 
present, see next Section) 

q L (x) = Z(f3)a A (P)q(x)+0(a 6 ) , (5) 

with a renormalization constant Z(f3) which is a finite function of the lattice bare cou- 
pling (3 = 2N c /gQ, approaching 1 as j3 — > oo. 

When defining the topological susceptibility, further renormalizations can appear. In- 
deed, already the continuum definition of Eq. © involves the product of two operators 
q(x) at the same point. Part of this contact term is necessary to make \ a positive quantity 
as requir ed by phenomeno logy because (q(x)q(O)) is negative for x ^ by reflection 

positivity l^ ' u l ' ' ~^ I Such a term is divergent, so that an appropriate prescription 
must be assigned to define it. This is easily done by making reference to Eq. ©, and it 
corresponds to fixing x = in the sector of zero topological charge. 
The lattice definition of the topological susceptibility 

XL = zZiq^qUO)) = ^ (6) 

X 

does not generally meet the continuum prescription, leading to the appearance of an 
additive renormalization M 



XL = Z{{5fa\{5) X + M(/3) 



(7) 



The quantity M{j3) contains the mixing with all local scalar operators appearing in the 
operator product expansion (OPE) of qL(x)qL{0) as x ~ in Eq. ©. 

The first lattice determinations - of x to °k account of the mixing with the identity 
operator (which gives the main contribution to M), but missed the multiplicative renor- 
malization, so that Z 2 x was measured, obtaining a value quite smaller than predicted 
by the Witten-Veneziano mechanism. Based on that, the idea was put forward that the 
field-theoretical discre tization of the topological charge might not be correct and the ge- 
ometric method ^ l ^ l , the cooling method ^^ l and Atiyah-Singer based methods'^ 
were developed. The field-theoretical method was then corrected by introducing Z and 

a correct subtraction M ^ESl The development of a non-perturbative tech nique, known 

as the heating method, for the numerical determination of these constants ^ ^ 1 finally 
brought about a reliable determination of x, free from the uncontrolled approximations 
involved in perturbation theory. 

The idea behind the heating method is that the UV fluctuations in qh{x), which are 
responsible for renormalizations, are effectively decoupled from the background topo- 
logical signal so that, starting from a classical configuration of fixed topological content, 
it is possible, by applying a few updating (heating) steps at the corresponding value of 0, 
to thermalize the UV fluctuations without altering the background topological content. 
This is surely true for high enough f3, i.e. approaching the continuum limit; in practice 
it turns out to be already true for the values of (3 usually chosen in Monte Carlo sim- 
ulations of gauge theories, being also favoured by the fact that topological modes have 
very large autocorrelation times, as compared to other non-topological observables (this 

autocorrelation time is particularly long in the case of full QCO^ai). 

One can thus create samples of configurations with a fixed topological content Q 
where the UV fluctuations are thermalized. Measurements of topological quantities on 
those samples can yield information about the renormalizations. For instance 

(Ql)q = Z{p) Q , (8) 

from which the value of Z(f3) can be inferred, while the expectation value of Q\ gives 

(Ql) Q = Z{f3) 2 Q 2 + V M{f3) , (9) 

where by V we intend now the four dimensional volume in lattice units and (-)q stands 
for the average within the given topological sector. 

To check that UV fluctuations have been thermalized, one looks for plateaux in quan- 
tities like (Ql)q or {Q\)q as a function of the heating steps performed: only configu- 
rations obtained after the plateau has been reached are included in the sample. Special 
care has to be paid to verify that during the heating procedure the background topological 
charge is left unchanged. This is usually done by performing a few cooling steps on a 
copy of the heated configuration and configurations where the background topological 

content has changed are discarded from the sample . Checks with other methods to 
measure the background topological charge (for instance based on operators that satisfy 
the Ginsparg-Wilson condition ^ | 27 | 28 ^ i ea( j to p er f ec tly consistent results^! 

A sample with Q = 1 can be used to measure Z and a sample with Q = (typically 
thermalized around the zero field configuration) can be used to determine M. Cross- 
checks can then be performed, using samples obtained starting from vario us confi gura- 

tions with the same or different values of Q, to test the validity of the method^^QI Once 



the renormalizations have been computed and the expectation value xl over the equilib- 
rium ensemble has been measured, the physical topological susceptibility \ is extracted 



= XL - m 
X a\f3)Z{(3) 2 K ' 

An analogous analysis and a similar method to compute the renormalization constants 
can be developed also for higher moments of the topological charge distribution!^. 

If renormalizations are large, i.e. if Z -C 1 and if M brings a good fraction of the 
whole signal xl, the determination of x via Eq. (ITOl can be affected by large statistical 
errors. A considerable improvement of the method is thus obtained by using operators 
for which the renormalization effects are reduced. This is the idea behind the definition 
of smeared operators^! q£ (x), which are constructed as in Eq. © but using, instead of 

the original links, the i-times smeared links ujp (x) defined as 

uW(x) = U^x), 

Ufix) = (l- c )U^(x) + C - £ Ut 1 Hx)Uj i i - 1 \x + a)Ut 1 \x + ^, 

Q = ±l 



Uf(x) 



±TxU\\x)\U\ r {x)\ 



where c is a free parameter which can be tuned to optimize the improvement. 



3 Topology at zero and finite temperature in QCD 

The method described in Section 2 has been used in the past in various applications; here 
we will briefly review the main results. The use of improved smeared operators has been 
essential to obtain high precision measurements, especially in the high temperature phase 
where the vanishing signal for x can g et completely lost in Eq. (flOl if renormalizations 
are large. Actually previous attempts to determine X across the transition have failed be- 
cause the errors involved in the determination were too large with respect to the physical 

signal. The second smearing level has revealed to be enough both in SU(2)&1 ( w i m 

c = 0.85) and SU(3p^(with c = 0.9) gauge theories. 

In Fig.^we display x at zero temperature in SU (3) pure gauge theory^las obtained 
from the measurements on 0-, 1- and 2-smeared operators at three different values of the 
inverse gauge coupling j3. The agreement among different operators (universality) and 
the good scaling to the continuum are apparent. From the 2-smeared operator the result 
X = (170(7) MeV) 4 is obtained if the value yfa = 420 MeV is considered together with 

the rations jy^F = 0.629(3). If instead the ratio T c /^/a = 0.646(7) is taken^then 
the result is x = (174(7) MeV) 4 . 

The behaviour of X across the finite temperature transition is an important ingredi- 
ent to understand the fate of the singlet axial symmetry at deconfinement and/or chiral 
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Figure 1: Topological susceptibility in SU(3) pure gauge theory calculated by using 1- and 2-smeared 
lattice topological charge operators (see text). The central band indicates the average result and its error bar. 
Left and right scales are shown in units of MeV and T c . 



SU(Nf) restoratioiPEl. The behaviour of x across the transition has been obtained^! 
and it is shown in Fig. |2] The susceptibility stays constant until the deconfinement tem- 
perature T c and then it undergoes an abrupt drop as shown in the Figure. Similar results 
are obtained in SU (2) pure gauge theoryE2 an d f or the unquenched theoryEH with 2 
and 4 staggered quarks. In Fig. |2 the behaviours of the quenched, Nf = 2 and Nf = 4 
theories with gauge group SU (3) are put together for comparison. In the unquenched 
calculation one has to take into account that the topological charge can also mix with 

other operators related to the anomaly^! Usually we neglect this mixing because it is 

rather small and surely smaller than the error bars from the simulation 1 ^. 

Another problem that has been studied by using the field-theoretical method is the 
possible spontaneous breaking of parity in Yang-Mills theory. An old theorem by Vafa 

and Witten precluded such a possibility 1 ^. However the authors implicitly assumed 
the breaking to be absent when they imposed that the derivative of the free energy with 
respect to 9 is zero at the minimum of the function — In Z(9). Our goal was to obtain 
a bound from lattice on the parity breaking effects, in particular to the electric dipole 
moment of neutral baryons: that was done by looking at a possible volume dependence 
of the topological susceptibility. After a thorough simulation on rather large volumes (up 
to 48 4 ) and huge statistics the following bound on the neutron electric dipole moment 

was obtained ^ (although it is still almost 5 orders of magnitude less precise than the 
corresponding experimental limit): d e < 3.5 10~ 21 e • cm . It must be stressed that the 
field-theoretical method enabled us to use such a large statistics and lattice sizes, a goal 
that would be hard to reproduce with other methods. 



4 Recent applications of the field-theoretical method 



In the applications of the Pisa method that we have reviewed in Section 3, a determination 
of x was usually needed at several values of the inverse gauge coupling f3, involving 
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Figure 2: Topological susceptibility in SU (3) pure gauge theory (left) calculated by using 1- and 2-smeared 
lattice topological charge operators (see text) at finite temperature across the transition at T c . The central 
band indicates the average result at zero temperature as obtained in Fig.Q A comparison is also made (right) 
with similar results obtained in full QCD with Nf = 2 and Nf — 4. 



several different determinations of the renormalization constants Z((5) and M((3). 

Suppose now that we want to study the behaviour of xi a ) on some parameter a, 
which does not affect the UV behaviour of the theory, while (3 and other bare parameters 
(like the masses of dynamical fermions, if present) are kept fixed. This happens in some 
interesting cases: we will consider the study of the behaviour of \ across the finite tem- 
perature deconfining transition on anisotropic lattices at fixed j3 and variable temporal 
extension (in this case a is the temporal extent N t of the lattice) and the determination of 
X across the finite density deconfining transition at fixed temperature (in this case a is the 
chemical potential /i). Renormalization constants are generated by quantum fluctuations 
at the UV scale, therefore they are independent of a and we can write: 

XL (a,P)-M(P) 

x{a) = ^mwf ■ ( } 

In the right hand side of Eq. (112b xl is the only quantity which depends on a: if one is not 
interested in an overall multiplicative factor, which is the case when studying the critical 
behaviour across a phase transition, the determination of a single additive renormalization 
constant is all that is needed to study the dependence of x on a - m those cases the 
field-theoretical method is much simpler and less computer-time consuming than other 
methods. 



4. 1 Topology on anisotropic lattices 

In the path integral approach to Quantum Field Theory, a finite temperature can be in- 
troduced by fixing a finite temporal extent r in the Euclidean space-time with (anti)- 
periodic boundary conditions for (fermionic) bosonic fields, the temperature being T = 
1/t. On a lattice this expression becomes T = 1/ (N t a{f3)). The temperature can thus 
be varied by changing either the number of lattice sites in the temporal direction or the 
lattice spacing, hence j3 (we are considering the pure gauge case, otherwise a would de- 
pend on the bare quark masses as well). Usually the second option is chosen: indeed for 




heating steps 

Figure 3: Comparison of the results obtained for (Q 2 l )q =0 /V as a function of the heating steps performed 
for the operator at the second smearing level on anisotropic lattices of two different temporal extensions: 

N t = 20 (circles) and Nt = 40 (triangles) . 

T near the QCD phase transition (T c ~ 200 MeV) and for the lattice spacings affordable 
by present computers (a^ 1 ~ 1 — 2 GeV) typical values of N t are too small to allow a 
fine tuning of T around T c by only varying the temporal length. 

The situation changes when using anisotropic lattices^! where different bare cou- 
plings are used in temporal and spatial planes, leading to different lattice spacings a s 
and at- One can thus use a very fine temporal spacing while leaving the spatial ones 
coarse, with an affordable computer time cost. Anisotropic lattices have been originally 
introduced to tackle problems related to heavy quarks, glueballs and high temperature 
thermodynamics which are not easily manageable otherwise. Here we will exploit the 
possibility of fine tuning T around T c by simply changing Nt and leaving a s and at 
fixed. 

We will consider the pure gauge plaquette action for N c = 3, 

S c = W~ £ ReTr (! " n ^)) + £ ReTr C 1 " n ^(*)) • d3) 

C ~ ^,i<j<3 c x ,i<3 

Both a s and at, as well as the renormalized anisotropy £ = a s /at, are functions of (3 and 
of the bare anisotropy 7. Several determinati ons of t hose parameters can be found in the 

literature, we will refer to the following value d p _ 6.25, 7 = 3.2552, which leads 
to £ = 4 and a s 0.021 fm. The critical temperature is reached for N t:C ~ 35, so that 
simulating at a different N t corresponds to T/T c = N tjC /N t . We have chosen a spatial 
lattice size N s = 24 which corresponds to a spatial extent of about 2 fm. 

We have determined x(T)/x{0) using the topological charge operator at the second 
smearing level and compared with previous results obtained on isotropic lattices. As 
a first step we have checked that the additive renormalization constant is indeed inde- 
pendent of the temporal extent Nt. We have computed M for two different temporal 
extensions, N t = 20 and N t = 40, corresponding respectively to T/T c w 1.75 and 
T/T c « 0.875. We have obtained M = (0.249 ± 0.010) • 10~ 5 for N t = 20 and 
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Figure 4: Topological susceptibility as a function of T/T c as obtained on anisotropic lattices by using the 
operator at the second smearing level (triangles) compared with previous results obtained using the same 

operator on isotropic lattice (circles). 

M = (0.242 ± 0.009) • 10~ 5 for N t = 40, i.e. a very good agreement which is also ap- 
parent directly from Fig. [3] where we report {Q\)q/V as a function of the heating step. 
An average of the two determinations of M has been used for other values of Nt. 

In Fig. H we report the quantity x(T)/x(0) = (xl(T) - M)/(xl(0) - M) obtained 
for the 2-smeared operator and compare it with the same quantity determined on isotropic 

lattices^! In the present case x(0) is not available and it has been fixed to the value 
obtained at the lowest available temperature, T = 0.875T C . A quite good agreement is 
visible. 

4.2 Topology at finite density 

An analogous situation is met when varying the chemical potential /i at constant temper- 
ature, which means at fixed (3, N t and quark mass values: indeed a finite fi does not affect 
the dynamics at the UV scale, at least until /x is small with respect to the UV cutoff (for 
very large values of \i Pauli blocking sets up with a consequent quenching of fermion 

dynamics at all scales). We report in this Section results obtained^ for the theory with 
iV c = 2, which is the only case where the sign problem does not make Monte Carlo sim- 
ulations at finite fj, unfeasible. A lattice of spatial size N s = 14 and temporal size N t = 6 
has been used, with dynamical staggered fermions corresponding to 8 continuum degen- 
erate flavours of bare mass am = 0.07 and an inverse coupling j3 = 1.5, corresponding 
to a temperature well below T c : a phase transition to deconfined matter is therefore ex- 
pected after increasing /j beyond some critical value \i c . The aim was to investigate the 
fate of topological excitations, hence of x> across the finite density phase transition. 

In the l.h.s. of Fig. [5] the estimates of M obtained for four different values of /i 
are shown. They look compatible within errors as they should be if the renormaliza- 
tion procedure is correct and no density effects are introduced into the subtractions. 
Hence also in this case a single renormalization constant is enough to determine the 
ratio x(/i)/x(0) = (xl(^) — M)/(xl{0) — M) which is displayed inther.h.s. ofFig.|5] 
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Figure 5: Additive renormalization (left) and x(/f)/x(0) ( rl ght) as a function of afi for the theory at finite 
density with 2 colours and 8 continuum dynamical flavours. 



It shows a clear drop at a critical value fj, c , analogously to what happens for the finite 
temperature phase transition. 

Acknowledgments 

The authors are greatly indebted to Adriano Di Giacomo for his acumen, the precious 
scientific collaboration and his sense of humanity. They also acknowledge collaboration 
with M.P. Lombardo and M. Pepe for the results obtained at finite density. 

References 

1. E. Witten, Nucl. Phys. B 156, 269 (1979). 

2. G. Veneziano, Nucl. Phys. B 159, 213 (1979). 

3. E. Shuryak, Comments Nucl. Part. Phys. 21, 235 (1994). 

4. M. Marchi and E. Meggiolaro, Nucl. Phys. B 665, 425 (2003). 

5. P. Costa et al., Phys. Rev. D 70, 116013 (2004). 

6. P. Di Vecchia, K. Fabricius, G. C. Rossi, and G. Veneziano, Nucl. Phys. B 192, 
392(1981). 

7. M. Campostrini, A. Di Giacomo and H. Panagopoulos, Phys. Lett. B 212, 206 
(1988). 

8. B. Alles and E. Vicari, Phys. Lett. B 268, 241 (1991). 

9. E. Seiler and I. O. Stamatescu, preprint MPI-PAE/PTh 10/87 (1987), unpublished. 

10. K. Osterwalder in "Constructive Field Theory", Lecture Notes in Physics, n° 25 
(1973), edited by G. Velo and A. S. Wightman. 

11. P. Menotti and A. Pelissetto, Nucl. Phys. Proc. Suppl. 4, 644 (1988). 

12. B. Alles, M. D'Elia, A. Di Giacomo and R. Kirchner, Phys. Rev. D 58, 114506 
(1998). 

13. I. Horvath et al., Phys. Lett. B 617, 49 (2005). 

14. E.-M. Ilgenfritz, K. Roller, Y. Roma, G. Schierholz, T. Streuer and V. Weinberg, 
hep-lat/05 12005 



15. M. Luscher, Comm. Math. Phys. 85, 39 (1982). 

16. P. Woit, Phys. Rev. Lett. 51, 638 (1983). 

17. E.-M. Ilgenfritz, M. L. Laursen, G. Schierholz, M. Miiller-Preussker and H. 
Schiller, Nucl. Phys. B 268, 693 (1986). 

18. M. Teper, Phys. Lett. B 171, 86 (1986). 

19. J. Smit and C. Vink, Nucl. Phys. B 286, 485 (1987). 

20. M. Campostrini, A. Di Giacomo H. Panagopoulos and E. Vicari, Nucl. Phys. B 
329, 683 (1990). 

21. A. Di Giacomo and E. Vicari, Phys. Lett. B 275, 429 (1992). 

22. B. Alles, M. Campostrini, A. Di Giacomo, Y. Gundiic and E. Vicari, Phys. Rev. D 
48,2284(1993). 

23. B. Alles, G. Boyd, M. D'Elia, A. Di Giacomo and E. Vicari, Phys. Lett. B 389, 
107 (1996). 

24. B. Alles et al., Phys. Rev. D 58, 071503 (1998). 

25. F. Farchioni and A. Papa, Nucl. Phys. B 431, 686 (1994). 

26. P. H. Ginsparg and K. G. Wilson, Phys. Rev. D 25, 2649 (1982). 

27. H. Neuberger, Phys. Lett. B 417, 141 (1998). 

28. M. Luscher, Phys. Lett. B 428, 342 (1998). 

29. B. Alles, M. D'Elia, A. Di Giacomo and C. Pica, Pro. Sci. LAT2005, 316 (2005), 
hep-lat/0509024 and work in progress. 

30. M. D'Elia, Nucl. Phys. B 661, 139 (2003). 

31. C. Christou, A. Di Giacomo, H. Panagopoulos and E. Vicari, Phys. Rev. D 53, 
2619 (1996). 

32. B. Alles, M. D'Elia and A. Di Giacomo, Phys. Lett. B 412, 119 (1997). 

33. B. Alles, M. D'Elia and A. Di Giacomo, Nucl. Phys. B 494, 281 (1997). 

34. G. Boyd et al., Nucl. Phys. B 469, 419 (1996). 

35. B. Lucini, M. Teper and U. Wenger, JHEP 0401:061 (2004). 

36. B. Alles, M. D'Elia and A. Di Giacomo, Phys. Lett. B 483, 139 (2000). 

37. D. Espriu and R. Tarrach, Z. Phys. C 16, 77 (1982). 

38. B. Alles, A. Di Giacomo, H. Panagopoulos and E. Vicari, Phys. Lett. B 350, 70 
(1995). 

39. C. Vafa and E. Witten, Phys. Rev. Lett. 53, 535 (1984). 

40. B. Alles, M. D'Elia and A. Di Giacomo, Phys. Rev. D 71, 034503 (2005). 

41. T.R. Klassen, Nucl. Phys. B 533, 557 (1998). 

42. N. Ishii, H. Suganuma and H. Matsufuru, Phys. Rev. D 66, 1 14507 (2002). 

43. B. Alles, M. D'Elia, M.P Lombardo and M. Pepe, Nucl. Phys. Proc. Suppl. 94, 
441 (2001); B. Alles, M. D'Elia and M.P. Lombardo, hep-lat/0602022 



